Two Critical velocities for a superfluid in a periodic potential 
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In contrast to a homogeneous superfluid which has only one critical velocity, there exist two critical 
velocities for a superfluid in a periodic potential. The first one, which we call inside critical velocity, 
is for a macroscopic impurity to move frictionlessly in the periodic superfluid system; the second, 
which is called trawler critical velocity, is the largest velocity of the lattice for the superfluidity to 
maintain. The result is relevant to the superfluidity observed in the Bose-Einstein condensate in an 
optical lattice and supersolid helium. 
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One characteristic feature of superfluidity is the exis- 
tence of critical velocity. According to Landau's theory of 
superfluidity this critical velocity is given by the speed 
of sound. In developing his theory Landau had his focus 
on a homogeneous superfluid, which has continued to be 
the focus of later studies 0, 0]. Experimental advances 
in recent years have brought to people's attention a dif- 
ferent type of systems, superfiuids in periodic potentials. 
They include a Bose-Einstein condensate (BEC) in an op- 
tical lattice Q and supersolid helium In the BEC sys- 
tem, the lattice is created by counter-propagating laser 
beams. For supersolid helium, the superfluid is defects 
(vacancies or interstitials) in the lattice self-assembled by 
helium atoms as commonly believed [H, @]- Interestingly, 
the crust of a neutron star can also be considered as a 
superfluid in a lattice [1, In addition, it is expected 
that the superfluidity of a paired Fermi gas in an optical 
lattice will be a subject of intensive investigation[lO(. 

In this Letter we show that the presence of the pe- 
riodic potential has non-trivial consequences, requiring 
a revisit of the concept of critical velocity. In contrast 
to the homogeneous superfluid which has only one crit- 
ical velocity, there are two distinct critical velocities for 
a superfluid in a periodic potential. The first one, which 
we call inside critical velocity, is for an impurity to move 
frictionlessly in the periodic superfluid system (Figflja)); 
the second, which is called trawler critical velocity, is 
the largest velocity of the lattice for the superfluidity to 
maintain (FigQJb)). These two critical velocities will be 
illustrated with a BEC in a one-dimensional optical lat- 
tice. 

The presence of the periodic potential plays a decisive 
role in the appearance of the two critical velocities. Be- 
cause of the addition of a periodic potential, two very 
different situations can arise in the superfluid system. 
The first situation is described in Fig[lja), where one 
macroscopic impurity moves inside the superfluid. The 
key feature in this situation is that there is no relative 
motion between the superfluid and the periodic potential. 
FigOJb) illustrates the other situation, where the lattice 
is slowly accelerated to a given velocity and there is a 
relative motion between the superfluid and the periodic 



potential. For these two different situations, naturally 
arise two critical velocities. 



^ superfluid 



^ superfluL 





FIG. 1: (a) A macroscopic obstacle moves with a velocity of v 
inside a superfluid residing in a periodic potential. The curly 
brace indicates that the superfluid and the periodic potential 
are "locked" together and there is no relative motion between 
them, (b) The lattice where a superfluid resides is slowly 
accelerated to a velocity of v. 



Both critical velocities can be measured with BECs in 
optical lattices. The inside critical velocity vi can be mea- 
sured with the same experimental setting as in Ref. Ill ] . 
where the superfluidity of a BEC was studied by moving 
a blue-detuned laser inside the BEC. For the trawler crit- 
ical velocity v t , one can repeat the experiment in Ref. fl2| 
where a BEC is loaded in a moving optical lattice. One 
only needs to shift his attention from dynamical insta- 
bility to superfluidity. For solid helium, only the trawler 
critical velocity can be measured. 

It is instructive to briefly review Landau's theory of 
critical velocity for superfluid before further discussion 
Consider a superfluid moving inside a small tube with a 
velocity of v as in FigO Suppose a single elementary 
excitation is generated and its energy is £o(p) and mo- 
mentum is p in the superfluid frame where the superfluid 
is at rest, then by the usual Galilean transformation, we 
have the excitation energy: 



e(v,p) =Eo(p) -v p. 



(1) 



in the reference frame where the tube is motionless. Lan- 
dau argues that if e(v, p) is negative for some excitations, 
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these excitations are energetically favored to be gener- 
ated and, therefore, the superfluidity is unstable. Since 
the low energy excitations are phonons £o(p) = U P, h is 
clear that e(v, p) can be negative only when v > u. This 
implies that the sound speed u is the critical velocity 
beyond which the superfluid becomes normal fluid. 
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FIG. 2: A superfluid moves inside a small tube with a velocity 
of v. 



all values of q and band index n, then the excitation of 
quasi-particles is not energetically favored and the sys- 
tem is a superfluid. Otherwise, the quasi-particles can 
be generated spontaneously and the liquid flow experi- 
ences viscosity. We can thus define a critical velocity, 
beyond which e' n (q) can be negative for some values of 
n, q. This critical velocity v% is given by 
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(6) 



We call Vi inside critical velocity to distinguish it from 
the other critical velocity that we shall discuss next. 



It is clear that the "tube" frame holds a special position 
in Landau's theory. It is only in this frame the tube does 
not disturb the system so that the normal fluid consisting 
of the excitation "gases" can be in a thermal equilibrium 
with the environment (l3j . However, in the presence of 
the periodic potential, the disturbance can be from either 
the external impurity (Fig. [Ha)) or the imperfection of 
the periodic potential (Fig.QJb)), and one needs to define 
the different "tube" frames for the different scenarios. 
For this reason and convenience of our future discussion, 
we call such a frame as thermodynamics frame. 

We focus first on the situation of FigQJa), where a 
macroscopic obstacle moves inside the periodic super- 
fluid with a velocity v. We follow closely Landau's 
argument and start our discussion in the superfluid 
frame where the superfluid is at rest. In such a peri- 
odic system, the elementary excitation (quasi-particle) 
is characterized by quasi-momentum q and the Bloch- 
band index n. If the excited state is described by the 
wavefunction , F„q, then its energy and momentum are, 
respectively, 



£„(q) = (*nq|-H"|*nq) ~ E , 
Pn(q) = (*»q|p|*nq), 



(2) 
(3) 



where H is the Hamiltonian of system with Eq being the 
ground state energy, p = Y]^ pj is the total momentum 
operator of the system. Note that without the periodic 
potential one would have simply p = Tiq. With the pe- 
riodic potential, this simple relation no longer holds. In 
the thermodynamic frame in which the obstacle is mo- 
tionless, the Hamiltonian of the system is transformed 
to: 



H' = H 



1 
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(4) 



where M is the total mass of the system. The excitation 
energy in this frame reads, 

4(v,q) = (*„q|i?'|*„ q )-^ =e„(q)-v-p„(q), (5) 

where E' = Eq + Mv 2 /2 is the ground state energy in 
the thermodynamics frame. This excitation energy de- 
termines the stability of the system. If it is positive for 
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FIG. 3: (a) Energy eo(q) and (b) momentum p(q) of a quasi- 
particle for a BEC in a one-dimensional optical lattice. The 
dashed line in (a) is for up(q) with u being the sound speed. 

We take a look at an example, a BEC in a one- 
dimensional optical lattice. We treat such a system with 
the standard Gross-Pitavskii equation, 



9 , 



1 d 2 



ijj + u(cosx)-0 + c\ip\ 2 ip . 



(7) 



In the above equation, the energy is in units of 
47r 2 ?i 2 m/a 2 , the unit of the momentum is 2?:h/a, and 
q is measured in units of it /a (a is the period of the op- 
tical lattice). For the details of the unit system, please 
refer to Ref. 14[. Following the procedure in Ref. 14[, 
we can compute its Bogoliubov excitations and, conse- 
quently, the energies and momenta of these excitations. 
Plotted in Fig|3] are the energy and momentum of the 
quasi-particle in this simple system. We see from this fig- 
ure that the dashed line for up(q) lies completely under 
e(q). This means that the critical velocity m is exactly 
the sound speed u in this simple example. 

The motion of an impurity inside a superfluid has 
been studied for a long time by moving ions inside liquid 
helium 15]. This technique was used to verify Landau's 
original criterion of superfluidity [lil] . which is impossible 
to verify by flowing superfluid helium inside a small tube 
as one would generate vortex and turbulence and destroy 
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superfluidity well before the superflow reaches the sound 
speed [3]. More recently, an obstacle created by a blue- 
detuned laser beam was moved back and forth inside a 
rather homogeneous BEC to test its superfluidity [TO] . We 
believe that this technique can also be used to measure 
the inside critical velocity for a BEC in an optical lattice. 
On the theoretical side, the motion of an impurity in a 
homogeneous superfluid was studied by Girardeau, who 
found the "onset of acoustical wave drag as the impurity 
speed reaches the speed of sound" [17 1 . 

We study now the case depicted in FigQJb), where the 
lattice is accelerated slowly to a given velocity v. In 
this case, the thermodynamics frame is the frame that 
moves along with the lattice. Viewing from the thermo- 
dynamics frame, the slow acceleration induces an adia- 
batic evolution of the Bose-Einstein condensation from 
the Bloch state at T point (k = 0) to a Bloch state with 
non- vanishing Bloch wave vector k = —mv/h [l8|, 
The same effect can be achieved by slowly turning on a 
moving lattice The excitation energy in the thermo- 
dynamics frame can be obtained in a similar way as that 
in Eq. ((2]), albeit the "ground state" for the present case 
has a macroscopic condensation at nonzero k. We thus 
have, 



(v,q) = (* nqik |.H'|tfnqik) - E (k) 



(8) 



where Eq(\l) is the energy of the ground state and l^nq^) 
denotes the excitation state with quasi-momentum q. We 
note that in the present case both the ground state and 
the excitation state have dependence on the condensa- 
tion momentum k = -mv/)i. As a result, the excitation 
energy e' n (v, q) depends on both the velocity v of the 
periodic potential and the quasi-momentum q of the ex- 
citation. The stability of the superfluid phase can be 
determined in the same spirit as that in the first case: 
if e^(v,q) can be negative for some finite q, then the 
superfluidity is lost. The critical velocity for this situ- 
ation is given by the smallest v such that e^(v,q) = 
for some finite q. We denote it as Vt and call it trawler 
critical velocity. This trawler critical velocity is related 
to the Landau instability discussed with the mean-field 
Gross-Pitaevskii equation in Ref.[2^, [2l|, (22|. Here in 
this Letter, this Landau instability is discussed in a more 
general setting. 

It is easy to demonstrate that the trawler critical ve- 
locity is different from the inside critical velocity. We 
consider a limiting case that the critical velocity is de- 
termined by the low energy phonon excitation (which 
means that Vt is necessary to be small). In this case, 
we can limit our consideration in the lowest Bloch band, 
and the full Hamiltonian of the system can be mapped 
to an effective one in which the boson moves in the free 
space (without the periodic potential) but with the renor- 
malized dispersion E(p) of the lowest Bloch band. For 
p ~ 0, -E'(p) ~ p 2 /2m* with m* being the effective mass. 



The effective Hamiltonian in the thermodynamics frame 
has the form 
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(9) 



where H is the effective Hamiltonian in the superfluid 
frame and H is the interaction potential projected in the 
lowest Bloch band. When an elementary excitation with 
quasi-momentum hq is generated, its energy is according 
to Eq.© 



e'(v, q) = £ (q) - h— v • q , 
m* 



(10) 



where £o(q) ~ uhq is the phonon excitation energy in 
the superfluid frame. It is clear that superfluidity is lost 
when v > m*uo/m. So, the critical velocity is 



-u , 



(11) 



which is different from the inside critical velocity Vi. Note 
that the velocity rrw/m* in Eq. (fTU|) is the group velocity 
of the superfluid. 
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FIG. 4: Inside critical velocity Vi and trawler critical velocity 
v t of a BEC in a one dimensional optical lattice as functions 
of the lattice strength v. (a) c = 0.02 (b) c = 0.1. The solid 
line is for Vi and the circles with dashed line represent Vt- 

We now demonstrate these two critical velocities, with 
the simple example, a BEC in a one-dimensional optical 
lattice. When this system is in the superfluid state, it 
can be well described by the mean-field Gross-Pitaevskii 
equation as in Eq.©. We can find both critical veloci- 
ties by numerically computing the Bogoliubov spectrum 
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of this system as in Ref. [LJ] . We have computed two dif- 
ferent cases: in one case the Bloch states corresponding 
to the trawler critical velocity are close to the T point; 
in the other these Bloch states are far away from the T 
point. For the first case, the results are shown in FigfJJa), 
where we see these two velocities have different trends: 
the inside critical velocities Vi decreases with the lattice 
strength while the trawler critical velocity increases. In 
the other case, both critical velocities decrease with the 
lattice strength. 

For supersolid helium, the critical velocity has not 
been measured up to date. The "critical velocity" (3.6- 
38/zm/s) measured by Kim and ( "n;m •">' is likely asso- 
ciated with the vortex quantization in the annular con- 
tainer. We are aware that the supersolid helium is a 
highly controversy topic at present [23j. 

In summary, we have tried to answer a simple ques- 
tion, "What is the critical velocity for a superfluid in a 
periodic potential?" The answer depends on the way how 
the critical velocity is probed. If you probe it by moving 
an impurity inside the BEC, you have the inside critical 
velocity; if you probe it by moving the lattice, you ob- 
tain the trawler critical velocity. These two velocities are 
different in nature and generally in values. For a BEC in 
an optical lattice, both critical velocities are measurable 
with current experimental techniques 11 , 12] • 

We thank A.L. Fetter drawing our attention to the 
early work on ion motion in liquid helium. We are sup- 
ported by the "BaiRen" program of the Chinese Academy 
of Sciences and the NSF of China (10504040,10604063) 
BW is also supported by the 973 project of China 
(2005CB724500,2006CB921400) . 



[2] L. D. Landau and E. M. Lifshitz, Statistical Physics (II) 
(Pergamon, Oxford, 1981). 

[3] P. Nozieres and D. Pines, The Theory of Quantum Liq- 
uids (II) (Addison- Wesley, New York, 1990). 

[4] O. Morsch and M. Oberthaler, Rev. Mod. Phys. 78, 179 
(2006). 

[5] E. Kim and M. H. W. Chan, Science 305, 1941 (2004). 
[6] A. F. Andreev and I. M. Lifshitz, Sov. Phys. JETP 29, 
1107 (1969). 

[7] A. J. Leggett, Phys. Rev. Lett. 25, 1543 (1970). 
[8] J. M. Lattimer and M. Prakash, Science 304, 536 (2004). 
[9] B. Carter, N. Chamel, and P. Haensel, Nucl. Phys. A 
759, 441 (2005). 
[10] M. W. Zwierlein, J. R. Abo-Shaeer, A. Schirotzek, C. H. 

Schunck, and W. Ketterle, Nature 435, 1047 (2005). 
[11] R. Onofrio, C. Raman, J. M. Vogels, J. R. Abo-Shaeer, 
A. P. Chikkatur, and W. Ketterle, Phys. Rev. Lett. 85, 
2228 (2000). 

[12] L. Fallani, L. D. Sarlo, J. E. Lye, M. Modugno, R. Saers, 
C. Fort, and M. Inguscio, Phys. Rev. Lett. 93, 140406 
(2004). 

[13] L. P. Pitaevskii, The Physics of Superconductors (I) 

(Springer- Verlag, 2003), chap. 2. 
[14] B. Wu and Q. Niu, New J. of Phys. 5, 104 (2003). 
[15] A. L. Fetter, The Physics of Liquid and Solid Helium (I) 

(Wiley, 1976), chap. 3. 
[16] F. Reif and L. Meyer, Phys. Rev. 119, 1164 (1960). 
[17] M. Girardeau, Phys. of Fluids 4, 279 (1961). 
[18] D.-I. Choi and Q. Niu, Phys. Rev. Lett. 82, 2022 (1999). 
[19] O. Morsch, J. Miiller, M. Cristiani, D. Ciampini, and 

E. Arimondo, Phys. Rev. Lett. 87, 140402 (2001). 
[20] B. Wu and Q. Niu, Phys. Rev. A 64, 061603 (2001). 
[21] A. Smerzi, A. Trombettoni, P. G. Kevrekidis, and A. R. 

Bishop, Phys. Rev. Lett. 89, 170402 (2002). 
[22] M. Modugno, C. Tozzo, and F. Dalfovo, Phys. Rev. A 

70, 043625 (2004). 
[23] D. Ceperley, Nature Phys. 2, 659 (2006). 



[1] L. D. Landau, J. Phys. U.S.S.R. 5, 71 (1941). 



